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In this paper, we establish two necessary conditions for a joint triangulation of two 
sets of n points in the plane and conjecture that they are sufficient. We show that these 
necessary conditions can be tested in 0(n 3 ) time. For the problem of a joint triangulation 
of two simple polygons of n vertices, we propose an 0(n 3 ) time algorithm for constructing 
a joint triangulation using dynamic programming. 

j_l ! 1 Introduction 



Let S be a set of points in the plane. A triangulation of S is a maximal set of line segments 
with endpoints in S such that no two segments intersect in their interior. A triangulation of 
S partitions the convex hull of S into regions not containing points in S that are bounded by 
triangles. Triangulating a set of unlabeled points in the plane under various constraints is a 
well studied problem in computational geometry [3, 4, 8]. 

Consider two sets A and B of points in the plane, where |A| = \B\ = n. Two triangulations 
T a of A and 7], of B are called joint triangulation (also called compatible triangulation) of A 
and B if there exists a bijection / between A and B such that (i) ijk is a triangle in T a if and 
only if f(i)f(j)f(k) is a triangle in and (ii) ijk and f(i)f(J)f(k) do not contain any point 
of A and B respectively (see Figure 1). The problem has applications in morphing [10, 11] and 
automated cartography [9]. 



1 The extended abstart of this paper appeared in the Proceedings of India- Taiwan Conference on Discrete 
Mathematics, Taipei, pp. 34-43, 2009. 



1 



(i) 
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(ii) 



Figure 1: Joint triangulations of two sets of points A and B: (i) bijection is not given, and (ii) 
bijection is given. 



The problem of joint triangulation of A and B has two variations depending upon whether 
the bijection between points of A and B are fixed in advance. The problem, where the bijection 
is not fixed in advance (see Figure l(i)), has been studied by Aichholzer et al. [1]. In this paper, 
we consider the other problem, where the bijection is fixed in advance (see Figure l(ii)). 

Let A = {ai, a 2 , . . . , a n } and B = {bi,b 2 , . . . , b n } be two disjoint sets of points in the plane, 
specified by their respective x and y coordinates. A line segment bfij is called the corresponding 
line segment of the line segment a^aj and vice versa. Similarly, a triangle bibjbk is called the 
corresponding triangle of aiCijCik and vice versa. Let T(A) and T{B) denote the set of all 
triangulations of A and B. The problem of joint triangulation of A and B, as stated earlier, is 
to find triangulations T(A) G T(A) and T(B) G T(B), if they exist, such that for each region 
bounded by a triangle aidjdk in T(A), the corresponding triangle bibjbk bounds a region in 
T(B) (see Figure l(ii)). The problem was posed in 1987 by Saalfeld [9], and since then, several 
researchers have worked on this problem but the problem is still open. 

The above definition of a joint triangulation of A and B needs some clarification. Consider 
triangulations T(A) and T(B) of point sets A and B respectively, shown in Figure 2. It can be 
seen that for every line segment aiOj in T(A), the corresponding line segment bibj is in T(B) 
and vice versa. However, the triangle (14(150,$ does not contain any point of A, whereas the 
corresponding triangle 6465^6 contains points of B. Thus the triangles bounding the regions 
are different and we do not consider this to be a joint triangulation. This gives rise to the 
definition of a component triangle as defined by Saalfeld [9]. A triangle in T(A) or T(B) is said 
to be a component triangle of the triangulation if it does not contain any point in its interior. 
Note that a triangle formed by three collinear points in A or B contains the middle point as 
its interior and therefore, such a triangle is not considered as a component triangle. Therefore, 
the problem of joint triangulation of A and B is to compute T(A) and T(B), if they exist, such 
that a triangle OiOjOk is a component triangle in T(A) if and only if the corresponding triangle 
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bibjb k is a component triangle in T(B). 




In the next section, we propose two necessary conditions for this problem and conjecture 
that they are sufficient. We also present an 0(n 3 ) time algorithm for testing these necessary 
conditions. If the given set of points A and B satisfy the two necessary conditions, we propose 
a greedy algorithm for constructing joint triangulations of A and B in Section 3. The proposed 
algorithm has been implemented and experimental results suggest that the algorithms correctly 
construct joint triangulations of A and B whenever A and B satisfy the two necessary con- 
ditions. Like two sets of points, a joint triangulation of two simple polygons of same number 
of vertices can be defined analogously. In Section 4, we present an 0(n 3 ) time algorithm for 
computing a joint triangulation of two simple polygons of n vertices. In Section 5, we conclude 
the paper with a few remarks. 

2 Necessary conditions 

Let CH(A) and CH(B) denote the boundary of convex hulls of A and B respectively. We state 
the first necessary condition for the existence of a joint triangulation of A and B, which relates 
the edges of CH(A) and CH(B), 

Necessary condition 1: If there exists a joint triangulation of A and B, then a^aj is an edge 
of CH(A) if and only if the corresponding edge bibj is an edge of CH(B). 

Proof: Assume on the contrary that there is a joint triangulation of A and B and an edge a^j 
is an edge in CH(A) but the corresponding edge bibj is not an edge in CH(B). Since ctiCij is an 
edge of CH(A), there exists only one component triangle (say, aia^a^) with (X% (X j clS £111 edge, in 
any triangulation of A. On the other hand, we know that any joint triangulation must include 
bibj in the triangulation of B. Since bibj is not an edge in CH(B) by assumption, there are 
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Figure 3: On the edge a 6 a 7 , a 6 a 7 a 8 and a e a 7 a 2 are successor triangles. The corresponding 
triangles b 6 b 7 b 8 and b 6 b 7 b 2 are also successor triangles on the edge b 6 b 7 



two component triangles (say, bibjbk and bibjbi) with bibj as an edge, in the triangulation of B. 
Since the component triangle OjOj-a/ is not present in the triangulation of A, this contradicts 
the definition of a joint triangulation. □ 

A triangle aiOjOk is said to be an empty triangle in A if it does not contain any point of A 
in its interior. Let Sa denote the set of all empty triangles in A whose corresponding triangles 
in B are empty triangles in B. Let Sb be the set of triangles corresponding to the triangles in 
Sa. It follows from the definition of a joint triangulation that only triangles from Sa and Sb 
can be component triangles in a joint triangulation of A and B. Let aidjCbk and ciiCLjai be two 
triangles in Sa such that they lie on opposite sides of their common edge ajOj. If bibjbk and 
bibjbi also lie on opposite sides of their common edge bibj, then aiCijdi is a called a successor 
triangle of a^ajau on the edge a^aj and vice versa. Analogously, bibjbi is also called a successor 
triangle of bibjbk on the edge bibj and vice versa. In Figure 3, a 6 a 7 a 8 and o 6 a 7 a 2 are successor 
triangles on the edge a§a 7 and their corresponding triangles b e b 7 bs and bsb 7 b 2 are also successor 
triangles on the edge b§b 7 . On the other hand, a%a 7 a% and a 2 a 7 a% are not successor triangles on 
the edge a 7 a% as a 2 and a$ lie on the same side of a 7 as- Since successors of aidjdk and bibjbk are 
defined jointly, in what follows, we say that ijl is a successor triangle of ijk on edge ij and vice 
versa. Observe that ijk can have more than one successor triangle on an edge ij. In Figure 3, 
(2,6,8), (7,6,8) and (3,6,8) are successor triangles of (5,6,8) on the edge (6,8). It is obvious 
that there is no successor triangle on any edge of the convex hull. 

Intuitively, if a triangle ijk is a component triangle in a joint triangulation, one of the 
successors on each edge of ijk that is not a convex hull edge is also a component triangle in the 
joint triangulation. Let S denote the maximal subset of triangles in Sa and Sb such that each 
triangle ijk in S has at least one successor triangle in S, on the edges ij, jk and ki that are 
not convex hull edges. Note that if a triangle ijk does not have a successor triangle on a non 
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convex hull edge, then ijk can not belong to S. We call triangles in S as legal triangles and S 
is called the set of legal triangles. Now, we state the second necessary condition. 

Necessary condition 2: If there exits a joint triangulation of A and B, then the set of legal 
triangles S is not empty. 

Proof: If there is a joint triangulation of A and B, then every component triangle in the joint 
triangulation has a successor triangle on each its non convex hull edges. So, every component 
triangle in a joint triangulation is a legal triangle and hence, the set of legal triangles S is not 
empty. □ 

Conjecture: There exists a joint triangulation of A and B if and only if A and B satisfy the 
two necessary conditions. 

Let us present an algorithm for testing the necessary conditions. The first necessary condi- 
tion can be tested by traversing the boundary of the convex hulls of A and B. Since the convex 
hulls can be computed in 0{n\ogn) time [3, 8], the first necessary condition can be tested in 
O(nlogn) time. 

For testing the second necessary condition, the algorithm starts by computing all empty 
triangles in A and B. It has been shown by Dobkin et. al [5] that all empty triangles in a set 
of n points in a plane can be computed in time proportional to the number of empty triangles 
which can be at most 0(n 3 ). So, Sa and Sb can be computed in 0(n 3 ) time. For every non- 
convex hull edge ij of all triangles in Sa and Sb, the algorithm checks whether there exists two 
triangles ijk and ijl on the edge ij in Sa as well as in Sb such that k and / lie on opposite 
sides of ij in both A and B. If ij satisfies this condition, then there are successor triangles on 
the edge ij. Otherwise, all triangles in Sa and Sb with ij as an edge are removed from Sa and 
Sb, and the remaining two edges of every deleted triangle are pushed into a queue Q. For each 
edge ef in Q, check whether there are successor triangles on ef. If the condition is satisfied, 
then ef is removed from the queue. Otherwise, all triangles in Sa and Sb with ef as an edge 
are removed from Sa and Sb, and the remaining edges of every deleted triangles are pushed 
into the queue Q. This process is repeated till either Sa and Sb become empty or the queue 
becomes empty. In the latter case, all remaining triangles in Sa and Sb have successors on all 
non-convex hull edges, in which case they form the set of legal triangles S. Note that that the 
cost of processing edges in Q can be assigned to deleted triangles which can be at most 0(n 3 ). 
We state the result in the following theorem. 

Theorem 1: Given two sets A and B of n points in the plane, the two necessary conditions 
for a joint triangulation of A and B can be tested in 0(n 3 ) time. 

3 An algorithm for constructing joint triangulations 

In this section, we present two algorithms for finding a joint triangulation of A and B which run 
in 0(n 3 ) time. We assume that the set of legal triangles S has been computed by the algorithm 
as mentioned in the previous section. If the set S is empty, clearly no joint triangulation exists. 
So, we consider the other case when S is not empty. 
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Figure 4: A joint triangulation of two simple polygons A and B. 



Constructing a joint triangulation of A and B involves finding a subset T of legal triangles 
in S forming a triangulation in A and the corresponding triangulation in B. The algorithm 
uses a greedy method to obtain T. Initialize S' = S and T = 0. Take any triangle ijk 
from S', add it to T and delete all triangles in S' that intersect the interior of the triangle 
ijk in either A or B. Repeat this process until S' becomes empty. Our claim is that the 
triangles in T form a joint triangulation of A and B. We have been unable to prove this claim, 
which would also prove the sufficiency of the two necessary conditions. On the other hand, we 
have observed experimentally that whenever S is not empty, the algorithm always finds a joint 
triangulation of A and B. Readers may use our software for experimentation, which is available 
at (http:/ / www.tcs.tifr. res. in/ ~ghosh/Joint-triangulation/joint-triangulation.html). 



4 Computing a joint triangulation of two simple polygons 

In this section, we present an 0(n 3 ) time algorithm for computing a joint triangulation of 
two simple polygons A = (a±, a 2 , . . . , a n ) and B — (b±, b 2 , ■ ■ ■ , b n ) using dynamic programming. 
Two points u and v in a simple polygon are said to be visible if the line segment uv lies totally 
inside the polygon. Let VG(A) denote the visibility graph of A, where vertices of A are vertices 
of VG(A) and two vertices in VG(A) are connected by an edge if and only if the corresponding 
vertices in A are visible in A [6]. The visibility graph VG(B) of B is defined analogously. We 
have the following observation (see Figure 4). 

Lemma 1: All edges of the triangles in a joint triangulation of A and B must belong to VG(A) 
and VG(B) respectively. 

Let IVG(A) denote the sub-graph of VG(A) such that an edge a^aj of VG(A) belongs to 
IVG(A) if and only if bibj is an edge of VG(B). Analogously, we define IVG(B). It follows 
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Figure 5: Testing the sub-polygon Qi,4 for a joint triangulation. 



from Lemma 1 that we have to consider only the edges of IVG(A) and IVG(B) in a joint 
triangulation of A and B. Since the visibility graph of a simple polygon can be computed in 
time proportional to the number of edges in the visibility graph, which can be at most 0(n 2 ) 
[7], IVG(A) and IVG(B) can be computed in 0(n 2 ) time. 

Let SUB(A) denote the set of all sub-polygons of A (including A itself) that can be formed 
by cutting A using only one diagonal of IVG(A). So, the size of sub-polygons in SUB(A) varies 
from 3 to n. We use a boolean function M{Q) to indicate whether a sub-polygon Q admits joint 
triangulation. Since all sub-polygons of three vertices in SUB (A) (say, <5i,3, Q2,3, ■ ■ •) admit 
joint triangulations as they are triangles, M(Qi^), M(Q 2^3), . . . are set to be true. Then the 
procedure considers sub-polygons <5i,4, Q2,4, ■ ■ ■ of SUB {A) having four vertices. Let Qi,4 = 
(dj, ai+i, ai+2, 0^+3) (see Figure 5). So, (lidi+z is the diagonal of IV G (A) used to cut A to 
form <5i,4- Let a k be a vertex of Qi,4 such that edges a^a^ and cika i+3 belong to IV G (A). If 
no such Vk exists, then set M(Q lt A to false. If a i+x = a k and the triangle (a i+i , a i+2 , a i+3 ) 
admits triangulation found in the previous step, then set M(Qi^) to true. If a i+2 = a k and 
the triangle (a*, a^+i, 0^+2) admits triangulation found in the previous step, then set M(Qi^) to 
true. Otherwise, set M(Qi^) to false. 

Similarly, the procedure considers sub-polygons Qi,5, ^2,5, • • • of SUB [A) having five vertices 
by locating all possible such vertices a k . This process is repeated till the sub-polygon of size n 
(i.e., A) is considered. In the following, we state the major steps of the procedure. 

Step 1: Divide A into sub-polygons using diagonals of IVG(A) to form SUB(A); 

Step 2: Consider each edge of A as a degenerated triangle; For each edge 0^+1 do 
M(aiai + i) := true; 

Step 3: For each sub-polygon Qj 3 of size three in SUB {A) do M(Qj j3 ) := true; size := 4; 
Step 4: For each sub-polygon Qj^ ize in SUB(A) do 

Step 4.1: If Qj t size = A then i :— 1, q := n, k := 2 and goto Step 4.3; 

Step 4.2: Let aja g be the diagonal used to cut A to form Qj tS i ze = (a,, a i+ i, . . . , a g ); 



fc := i + 1; 
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Step 4.3: If ' a^ak and a q cik are edges in IVG(A) and two sub-polygons formed by removing 
the triangle (ai,a,k,a q ) from Qj tS i Ze admit joint triangulations then 
M(Q jtSize ) := true; 

Step 4.4: If k ^ q — 1 then k := k + 1 and goto Step 4.3; 

Step 5: If size ^ n then size := size + 1 and goto Step 4; 

Step 6: If M(A) is true then by backtracking identify diagonals of IVG(A) giving a joint 
triangulation else report that there is no joint triangulation. 

Step 7: Stop. 

Since the procedure uses triangles formed by edges of IVG(A) and IVG(B), and these 
triangles are added one at a time (i.e., aja^a^) to verify whether a joint triangulation exists 
for the sub-polygons formed by the union of triangles verified so far, the procedure correctly 
computes a joint triangulation of A and B if it exists. Since the number of sub-polygons in 
SUB (A) can be at most 0{n 2 ) and the procedure can take 0(n) time for testing each sub- 
polygon, the overall time required by the algorithm is 0(n 3 ). We state the result in the following 
theorem. 

Theorem 2: Given two simple polygons A and B of n points, a joint triangulation of A and 
B can be constructed in 0(n 3 ) time. 

5 Concluding remarks 

Let us mention some extensions of the basic problem. An immediate extension is to find a joint 
triangulation of k sets of labeled points. It is easy to verify that for such a joint triangulation to 
exist, boundary of the convex hulls of all sets of points must contain the same edges. Further, 
the notion of a successor triangle can be extended to any number of sets of points in a natural 
way. A triangle ijl is a successor of a triangle ijk on the edge ij if and only if it is a successor 
in all point sets. Thus we may define the set of legal triangles in an analogous way. We believe 
that the same conjecture holds for any number of sets of points. 

Further generalizations are possible by considering triangulations of objects other than just 
point sets. In particular, we can consider triangulations of any connected polygonal region with 
points and polygonal holes inside. The only difference here is that a triangle containing an edge 
of a hole boundary may not have a successor on that edge. Thus one necessary condition is 
that the hole boundaries must contain the same set of edges in all point sets. The definition of 
a successor triangle and a legal triangle may be modified accordingly, and the same algorithms 
can also be used. Again, we have observed empirically that if the set of legal triangles is not 
empty, there exists a joint triangulation, and it may be constructed in the same greedy fashion 
as for two point sets. 

If there is no joint triangulation of A and B, it may still be possible to obtain a joint 
triangulation by adding some points (say, m Steiner points) in A and B. Naturally, it is 
desireable to add the smallest m so that A and B admit joint triangulation. Aichholzer et 
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al. [1] showed that it is always possible to obtain joint triangulation of A and B (without 
a bijection) by adding a linear number of Steiner points. One would expect a better bound 
where bijection between A and B is given in advance. In the case of simple polygons A and B 
(without a bijection), Aronov et al. [2] showed that an addition of quadratic number of Steiner 
points is sufficient and sometime necessary for constructing a joint triangulation. 
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